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I INIR,Q.pyCTION 

This report describes the work performed during the second semi- 
annual period under contract NGL-39-01 0-001 . The work performed 
is partially a continuation and completion of work started in the first 
reporting period. Besides, new work has been started primarily in the 
direction of gaining additional understanding of the theoretical relation- 
ships determining the current-voltage characteristic of silicon solar 
cells . 

The calculations performed in the first reporting period to gen- 
erate design data for the light generated current in dependence on 
material and physical device parameters have lead to the statement in 
the report covering the first period, that the data might lend themselves 
to the preparation of a "solar cell design handbook" . Although the de- 
sign data obtained during that period fulfilled the requirements of the 
present contract, it was felt desireable to rework these data and aug- 
ment them by additional data to form the beginning of such a design 
handbook. Also, it was thought desireable to appropriately explain the 
mathematical basis on which such design curves are based. This has 
been done, and section 2 of this report contains essentially the first 
section of this "design handbook", dealing with the collection effi- 
ciency part of the design tasks. 

Since certain recommendations arose from the design calculations 
carried out relative to the collection efficiency, it was found also 
desireable to first investigate the effect of such structure changes 
on the current-voltage characteristic. This investigation has been 
carried out, and the existing theoretical relationships for diffusion 
current contributions from the base region and the diffused region have 
been expanded to cover the two-layer model including drift fields , 
analogous to that used for the calculations of collection efficiency. 
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A new computer program has been generated based on 
the new theoretical relationships , and design calculations have 
been carried out. Although the evaluations of these calculations 
is not complete , preliminary investigation of the data indicates 
that any structural and material parameter measures taken to im- 
prove the collection efficiency also reduce the dark diode current, 
and thus improve the current- voltage characteristic. 

A computer program for the reduction of experimental data on 
the current- voltage characteristic of solar cells has been described 
in the first semi-annual report. Further work on this program has 
been carried out, and has led to successful completion. The pro- 
gram has been found extremely useful, saving not only considerable 
amounts of time in the reduction of such data , but also providing an 
accuracy in determining the 5 constants of the two-exponential current- 
voltage characteristic, not possible by any other method. Use of this 
program on a relatively small number of cells so far has led to the 
conclusion that the two- exponential relationship appears to be the 
appropriate one to describe the current- voltage characteristic of 
silicon solar cells . Comments on the practice of reduction of computer 
data on the experimental current- voltage characteristic, a description 
of the computer program and of its use, and some results obtained with 
it are discussed in section 4 of this report. 

The experimental effort towards verification of the collection efficiency 
improvement in the short wavelength region by reduction of the surface 
recombination velocity and reduction of the diffused region thickness has 
been delayed, since a similar program has been carried out at Comsat 
Laboratories. The program there has been extremely successful, and has 
led to verification of the predictions made in the beginnings of this program. 
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II 


THE "SILICON^ SOLAR'' CELL DESIGN’* H A Nl} BOOK" 


1 Introduction 

This is the first part of an endeavor to create a "silicon solar cell design 
handbook". Through such a handbook, it should become possible to reduce 
solar cell design from an art and science to an engineering task. This objec- 
tive can be accomplished by either of two methods: in the old-fashioned way 
of presenting design curves and data tables, and the modem way of computerized 
design. Here, a time-inverted approach is used, primarily for instructional pur- 
poses: the computer has been used to provide the data from which design curves 
have been prepared. From the design curves, it is easier to gleen the influence 
of the various design parameters, than it would be from the results of a computer 
design program which yields the optimized design data directly i Getting a feel 

for the influence of the various material and dimensional parameters of the device, 
and of their interrelationships is still important at the present stage of develop- 
ment. 

A design process in general consists of an analysis and a synthesis. In the 
analysis part, the influence of material and design parameters on each component, 
subsystem, or performance attribute is established. In the synthesis part, these 
elements are combined to provide the overall performance of the system, possibly 
necessitating recycling through the first stage to provide an optimized design. 

The design approach is logically split into three parts, concerned with the 
light generated current , the open circuit voltage, and the fill factor. A principal 
measure of performance relative to the light generated current is the collection effi- 
ciency. The collection efficiency is strongly influenced by the spectral distribution 
of the incoming light*. The following design data for collection efficiency are re- 
stricted to airmass zero sunlight. Since open circuit voltage and fill factor are not 
dependent on the spectral distribution, at least in the commonly used silicon solar 
cells, the design data relating to these two quantities are of broader validity. 
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2 . The Basis for the Collection Efficiency Design Curves 


In order to perform such an analysis and synthesis, it appears most practical 
to break the total task into logical parts. This can be done by starting from the 
governing performance criterion, the efficiency for solar energy conversion, fol- 
lowing well established practices for the most part. Overall conversion efficiency 
of a solar cell, like any efficiency, is determined by the relationship: 


where: 
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integrated over all wavelengths comprising the intensity spectrum of the light inci- 
dent upon the cell, (ref.l) The maximum output power available is the product d£ 
a current density and a voltage: 


P . = j .V 
out max max 


( 3 ) 


Here , all power values and the current density have been referred to unit area of 
light-exposed cell-surface. The power output has been related to the light gener- 
ated current and the energy gap by means of two factors: 
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( V.F.), called the "voltage factor 1 ’, and ( F.F.), called the "fill factor", are 
largely determined by the characteristic of the rectifying potential barrier of the 
device ( in the Si solar cell*, a pn-junction). However, shunt current and series 
resistance effects can degrade both factors below those which are based on the 
characteristics of the barrier itself. As a result of this observation, the voltage 


factor is best described as a product: (V.F. ) = ( V.F.) 


J 


(V.F. ) 


A. 
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The subscript "J" refers to the voltage factor as it would be obtained in an "ideal" 
solar cell.,' where only the characteristic of the rectifying potential barrier of the 
device determines the open circuit voltage. However, artefacts such as process 
or ddsign effects which result in effective shunt currents, can lower the open cir- 
cuit voltage. These effects are lumped in the second voltage factor designated 
by subscript "A" . In silicon solar cells , the shunt currents have generally been 
reduced to the point where they -have negligible influence on the open circuit volt- 
age. Therefore, (V.F. ) <&l, and only the first factor needs consideration. In 
this case, the subscript "J" is generally omitted, and the voltage factor understood 
to be based on the junction characteristic only. 


The Fill- Factor ( F.F. ) can similarly be broken into 2 factors: 

(F.F.) = (C.F.) ' (A.F.) (4c) 

where the “Curve Factor" (C.F) has been used in the literature (ref.l) to designate 
the part based on the characteristics of the potential barrier, and where (A.F.) is 
introduced herewith to account for additional "softening" of the IV - characteristic, 


as it is caused by artefacts such as shunting or series resistance. 

1 1 = - 


Thus; 


E q F (J) G (A) 


(5) 


m 


where: F (J) = (V.F.) ' (C.F.) is a function dominated by the potential barrier, 
being less strongly dependent on j^, and G (A) = (V.F.)^ . (A.F.) is a function 
primarily of the artefacts, while J L depends primarily on the material properties 
of the diffused and the base regions of the device, and rather little on the para- 
meters of the transition region associated with the potential barrier. 


It is: 


j L = q 
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where ^s(A) is called the overall collection efficiency, since it relates the 
number of charge carriers collected and made utilizable in current flow through 
on- outside circuit to the number of photons available for conversion. 
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It is: 


or: 


P ln (X) d\= N ph (X ) 
a N 


E ph (X ) A ; 
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where N , ( \ ) is the number of photons incident in the wavelength range dA 

ph 

around wavelength X, h being Planck's constant and c the velocity of light. 

E ( A ) is the energy of the individual photon of wavelength A. 
ph 
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Not all of the incident energy is absorbed in the cell: some is reflected 
from the surface, another part is not absorbed, but effectively transmitted 
through the device. The reflected energy is: 

- r(X) : ( 9 ) 

where r(X ) is the reflection coefficient at wavelength A. The transmitted 
energy is given by Lambert's law of absorption: 
r, r \ ^ \ d; 

P^A) = (A) e 


(10) 


where d is the thickness of the absorbing layer, P q the energy entering the layer, 
and<x(A) the absorption coefficient at wavelength A* for simplicity, it is as- 
sumed that reflection at the backsurface does not take place. 

Them P in (X ) [l - r (X )] [ 1 - e" * (X) ^JdX (11) 

is the energy actually absorbed in the layer of thickness d in the wavelength in- 
terval dA at A. Both relations - for reflection and for absorption are generally 

used for light intensities , or incident energy fluxP . It is seen that they can also 
be applied to N 
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Here, 1| c ( A ) is £ he collection efficiency of the device, related to the number 

of photons actually absorbed within it. Since useful absorption does not take place 
beyond the absorption edge, that is absorption with generation of minority carriers, 
it suffices to extend the integral at its upper limit to a cutoff wavelength A^. deter- 
mined by the energy gap of the semiconductor used: 

A = — ; (13) 

G ' 


The expression for the total light generated current j can thus be written in the 

L 


j L (A )dA; 



At this point the first assumption will be introduced. It is possible that 

(A ) is not only a function of A, but also of N , (A*)/ where A* can be any 

ph 


wavelength other than A within or outside of the range of integration. If this is 


the case, the integral is not linear m N (A). This will be the case, for in- 

ph 

stance, if trapping of charge carriers occurs to a significant degree , or if the 


total incident light intensity is so large as to generate a minionty carrier 
density m any region of the solar cell which is not small compared to the major- 
ity carrier density, so that minority carrier lifetime can no longer be considered 


independent of minority carrier density and, m general, the normal assumptions 
of small signal theory can not be made. Fortunately, in normal silicon solar 
cells, and probably m many other types of solar cells, the conditions for valid- 


ity of the assumption are fulfilled, even for the use of rather large concentration 
ratios. A device which is known to violate the assumption is the Cu S-CdS solar 

Li 

cell. It will therefore be necessary to test the validity of the assumption before 
applying the relationships to be subsequently developed herein to other types of 
solar cells . Although the general methods applied here may be useful in these 
other cases also, modifications of some details will be necessary. 

It is practical at this point to introduce weighted average values for the 
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reflected energy, the energy actually absorbed m the wafer, and the collection 
efficiency 7^ 



The fraction of the energy absorbed m the device to that entering it is given by 



while the reflected energy is: 



Since the reflection can, in practice, be reduced to very small values, R ^ 
is assumed to be zero in the following. However, actual reflection characteris- 
tics can very readily be introduced into the subsequent equations by replacing 

Nph (A ) with N ^ (A) [l -r(X)J . Note; if 1 - r (A ) is significant and 

has a strong dependence on wavelength, the design curves given in the following 
will not be accurate and should be re-computed for the specific reflection charac- 
teristic. 
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3. Collection from the 3 major regions of the ce ll 


Since contributions to the light generated current j can be identified as 

Li 

originating m three different regions of the solar cell, it is practical to write: 


J L = / [ j L,D <A> +i L,T ?W+) L,B SA, J dA 1 


(19) 


where the subscripts D,T, and B refer to the diffused region, transition region, 

and base region of the cell, respectively. In the following, it is tacitly assumed, 

without loss of generality, that the diffused region is the one adjacent to the light 

exposed surface. If an "inverted" structure should be considered, it is merely 

necessary to enter the parameters of the diffused region intol^^^^ ^ and those 

of the nohdiffused region into 71 ,, ^ 

^coll,D. 


Three collection efficiencies can then be introduced, one each connected with 
one of the major regions. 

i-« N Ph a>-j r [i-e-«a > x TF j 

0 + [e-*< A > *T, F *T. R J 

+ [e-* (X) -e-* (A)d ] ^coII.e}^ * (20 > 

Here', x^ ^ is the location of the front boundary of the transition region 
(depletion region), that is the one situated more closely to the light exposed 
surface of the device, to which x = 0 has been assigned. A second assumption 
has been made here, containing the condition that the entire structure is plane- 
parallel so that all surfaces of significance are planes and are parallel to the 
light exposed surface. Not included m this parallelity condition are obviously 
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the surfaces which form the edges of the wafer. 

x is the location of the rear boundary surface of the transition region. 

T,R 

Equation (20) can be written as three separate integrals, designating the light 
generated current contributions from the three regions: 



and: 



fl/ G N oh f 1 - e ~* a> X T.rJ dA 
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Similarly: 
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and 
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With this, now /the foundation has been laid for the presentation of design 
data for the various aspects or parts of the device. R permits evaluation of 

cl.DS 

the direct effect of wafer thickness d on the overall collection efficiency. Intro- 
duction of R d , R g , and R T permits separate optimization of the mam regions of 
the solar cell, the design of which can be carried out relatively independently. 

The performance criteria here are the collection efficiencies for the three regions. 
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4 . The Collection Efficiency Design Curves 


Basic to the use of the collection efficiency design curves is the knowledge 

of R , , R^, R m/ and R_ . R , is given in Fig.II-I as function of wafer thickness 
abs D T B abs 

d. R and R are directly obtained from Fig.II-2, when the transition region thick- 
D B 

ness of the pn-junction is taken to be zero. If the transition region thickness is not 

zero, then R and R are also obtained from Fig.II-2, but in this case, by using x. p 
D B J r* 

for reading off R^, instead of x^ as the parameter,, and similarly Xj ^ instead of for 

obtaining R . 

B 


R^ is then obtained from 


R T 1 " r d r b ; 


(24) 


The collection efficiency of the transition region has been taken as 100% up 
to this point. 

It is -important for the user of the Design Handbook to realize the restrictions 
and .as sumptions , under which the data have been prepared. It has to be left to him 
to assure himself, before using the data, that his case at hand does not significantly 
violate these assumptions. This is the reason why, in the discussion of the previous 
sub-section, the assumptions made have been stressed, and why, in the following 
paragraphs, further discussions follow on the assumptions made and on their impact 
on the design data. 

The weighted average values of the collection efficiency and the ratio 
absorbed, as represented by equations (22a)-(24a) and (22b)-(24b), respectively, 
would ideally be obtained from an integration process as shown ,m these equations. 
In actuality, the integration has been replaced by an addition process over 50 nm 
wide wavelength intervals, as outlined in section II of the first semi-annual report. 
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a single pass through a silicon layer of thickness 
d to the number entering the layer, as contained 
in alrmass zero sunlight up to 1.125 jam wavelength. 





Wafer Thickness d (cm) 


H-2 Ratios Rj-j and R g of the number of photons of airmaas zero 

sunlight which are absorbed In the diffused region and the 
base region, respectively, as function of wafer thickness 
d. The ratio Is taken relative to the number of photons 
absorbed In the wafer of thickness d, and Is presented 
parametlcally for various thicknesses x, of the diffused 
region (depletion region width zero). -Where the depletion 
region width Is not zero, x_ p and x ^ will be used for the 
diffused region and the base' region, \tespectively. 




The determination of the weighted average collection efficiency requires 
a knowledge of the collection efficiency as function of wavelength, that is^A). 

These functions ?£(A ) have been calculated by use of -the computer program 
SOCEL2, which was discussed m section II of the first semi-annual report. This 
program is based on an exact solution of the continuity equation for minority carri- 
ers . The only basic assumption entering into this differential equation is related 
to the recombination statistics and is contained in the concept of minority carrier 
lifetime to be independent of the concentration of the carriers or the way in which 
they are produced. This is related to the fundamental assumption of small signal 
theory, that conditions never deviate substantially from thermal equilibrium con- 
ditions. This entails, for instance, that the injected carrier density never sig- 
nificantly changes the majority carrier concentration. 

An additional, more restrictive assumption has been made to obtain analytical 
solutions from the continuity equation. This assumption requires the independence 
of the material constants from the space coordinate x. In those cases where no im- 
purity density gradient exists within a region, this assumption is fulfilled. However, 
in general, it is not fulfilled in the diffused region. It is also not fulfilled m any 
subregion where a drift field is purposely introduced. Here, the minority carrier life- 
time % and mobility will be functions of x. The drift field E depends on the impuri- 
ty density distribution with x. E is independent of x, if this impurity density 
distribution has an exponential dependence on x. 

The method used by this author was to assume constant material parameters 
V and yd, approximately equal to the weighted average value for the region. Thus, 
the larger part of the region has material constants of slightly pessimistic value, 
while a smaller part uses a slightly optimistic value. 

As mentioned, constancy of the drift field E requires either a constant or an 
exponential impurity distribution. Complementary error function distributions are 
more likely to be encountered. However, an approximation of this distribution by 
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an exponential function is reasonably close for the major part of the distribution, 
and thus the assumption of constant drift field is frequently not a bad one. 

This author has found earlier (ref. 2) that the same performance can be ob- 
tained from the diffused region either by minority carrier diffusion only, with a 
certain diffusion length, or by drift and diffusion combined, with a shorter diffusion 
length and an electric field corresponding to the expected impurity diffusion profile. 
Since neither electric field nor diffusion length are measured diractly, but are in- 
ferred from the performance, this author has preferred to use the simpler case of 
no field and longer diffusion length. Once more is known about the actual diffusion 
profiles, the design curves should be reviewed. In the meantime, they form good 
approximations . 

In the same paper (ref. 2), this author used the same method used here to 
investigate the influence of a drift field m the base region on cell performance 
under the influence of damaging nuclear particle- radiation. The- drift field was used 
m a subregion near the pn-junction. Several other authors have since attempted to 
improve on this method by analyzing the same problem. The authors of ref. 3 have 
included a linear dependence of /A, X , and E on x in the first subregion of the 
base. However, they considered the second subregion in the base as a passive 
layer, that is, not contributing any carriers to the collection process. The effect 
of this assumption is to move the optimum to greater subregion 1 thickness, and to 
make the analysis not equivalent to that of ref. 2. They further assumed the base 
region to be infinitely thick and thus could neglect any effect of surface recombi- 
nation velocity on the back surface. 

The authors of ref. 4 went to much greater pain in utilizing the known de- 
pendence of lifetime and mobility on impurity concentration, and in using an electric 
field distribution corresponding to the complementary error function distribution of 
the impurities . Also, they do not appear to have made the assumption of non-col- 
lection from the subregion behind the drift field region. These authors found that 
there is "only a slight reduction of collection efficiency in going from the expo- 
nential case to the complementary error function case, if the width of the drift field 
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in the exponential case is taken to be equal to the thickness of the epitaxial lay^" 
which essentially means to exclude the high concentration portion of the comple- 
mentary error function distribution. In general, the results of ref. 4 appear to be 
extremely close to those obtained in ref. 2 by the "constant field and material 
parameter" assumption. More recently, another pair of authors (ref. 5) has dealt 
with the same subject, primarily trying to reconcile the differences between refer- 
ences 3 and 4. However, their approach is also questionable since they neglected 
generation of electrons in the base subregion behind the drift field region. Further, 
use of the boundary condition of continuity of n and ~ at the subregion interface 
expresses the physical requirement of continuity of minority carrier flow across this 
interface only, if the drift field at this interface is zero on both sides of the inter- 
face. However, ref. 5 treats cases where this condition is not fulfilled . Here, the 
boundary conditions of ref. 2 would have been required. Despite these discrepancies, 
the results of their computations indicate, that the assumptions of constant material 
parameters and constant drift field do not substantially alter the dependence of short 
circuit current or collection efficiency on the physical parameters studi»d. More- 
over, they show that the constancy assumption even leads to the proper quantitative 
values, if judicious choice of the constant values is made. Analysis of the data 

shows that despite all the differences in assumptions between ref. 2 to ref. 5, the 
quantitive results do not appear to differ by more than 3 to 5 percent in comparable 
cases, and thes,e differences. are probably-at leSs impartially caused by differences 
in the mobility versus impurity concentration data used. 

Another question occasionally raised bohcerns'the integration over wave- 
length. Ref. 3 to ref. 5 have followed Kleinman's approach (ref. 6) of integrating 
first over wavelength, thus obtaining a non-analytical forcing function in dependence 
on x for entry into the ordinary, linear, inhomogeneous differential equation. This 
author solved the differential equation analytically with constant coefficients for 
discrete wavelength values in the range of interest, and added the solutions to ob- 
tain the total light generated current. The latter method has the advantage that 
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j (X) and ^ A) are obtained, that is the spectral response. This permits more 
detailed comparison with experiment. While analytical solutions are not available 
with non-constant coefficients, the numerical methods could still be applied to 
obtain spectral response data. 

Both methods , that is integration over wavelength either before or after 
solution of the differential equation, are mathematically equally valid. In either 
case, the integration, due to the non-analytical nature of N ^ (A)/ has to be re- 
placed by a summation with adequately fine steps . Thus , the forcing function is 
a sum. Since we deal with an ordinary, linear differential equation, the solution 
of the inhomogeneous equation is the sum of the general solution of the homogeneous 
equation and a particular solution of the inhomogeneous equation. Since the forcing 
function is a sum, the particular solution of the inhomogeneous equation can be ex- 
pressed as a sum of the particular solutions for each of the addends of the forcing 
function. This is the principle of the “spectral response" approach. Both approaches, 
if properly executed should yield the same results. 

To summarize the discussion on the limitations of the data used: 

1. The data relating to homogeneous subregions are accurate. 

2. Where subregions with drift fields have been analyzed, an approximation has 
been used. From results obtained by other methods, it may be concluded 
that: 

a . The relationships between various design parameters and performance 
have been obtained in very close approximation. 

b. Absolute values of collection efficiency may be off by a few percent, 
depending primarily on the choice of material parameters , such as 
mobility, made. 

After this preamble, then, the design curves for collection efficiency from 
the diffused region and the base region are presented m the following . The meaning 
and impact of various design approaches have been discussed in the first semi-annual 
report, and the curves are rather self-explanatory, so that further discussions at this 
point should not be necessary. The data presented previously have been augmented 


- 19 - 



by new collection efficiency curves for the base region for diffused region thickness 
of 0.8, 2.5, and 6.4 yUm, to complete the set of design parameters covered and to 
round out the visibility of impact of variation of the individual parameters. With in- 
creasing thickness of the diffused region, the light generated current contribution 
from the base region decreases rapidly, but the collection efficiency from both the 
diffused and the base regions decrease also. 
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Diffused Region Collection Efficiency 7 C oll,D 


x T>F = 02/i m 


s= I0 2 cm s' 



Fig.n - 3 Collection Efficiency ^ from the diffused region as function 

of minority carrier lifetime , for several values of surface recom- 
bination velocity s. (Diffused region thickness Xj p = 0.2/rm.) 



Diffused Region Collection Efficiency 

7coll, D 


'*t,f = 0.2/zm 



Fig. U-4 Collection Efficiency D from the diffused region as function 

of surface recombination velocity, for several values of minority carrier 
lifetime*^. (Diffused region thickness Jtj p = 0.2^m.) 


Diffused Region Collection Efficiency 7coll|D 



Fig. n - S Collection efficlency^,^ D from the diffused region as 

function of minority carrier lifetime , for several values of 
surface recombination velocity s, (Diffused region thickness 
* T F - 0.4/fm). 
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Fig. IX - 6 Collection efflclency^^ Q from the diffused region as 
function of surface recombination velocity, for several 
values of minority carrier lifetime?^ . (Diffused region 
thickness^ p = 0.4ytm] 
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Fig. 11-9 Collection efficlency^^ ^ from the diffused region as function 

oE surface recombination velocity, for several values of diffused 

region thickness _ with minority carrier lifetime values T 
I * ‘ P 

which yield maximum light generated current. 
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Fig. n - 10 Design curves for the material parameters as function of the 
diffused region thickness. The curves yield the values of 
minority carrier lifetime and surface recombination velocity s , 
for which a collection efaciency of 97.5%, 95%, and 90%, 
respectively, can be obtained, assuming the other one of these 
’ 2 parameters being adjusted so that this Is possible. 
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Fig. n- 14 Collection Efficiency 7J coll B from the base region as function 
of wafer thickness d, for 3 values of diffusion length In the base 
region, with and without a drift field in front of the back contact. 
Airmasa zero sunlight, T = 300° K. Xj R = 0.8/«m. 
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Fig, 11-19 Collection efHciencj^^^j from the base region as function 

v of thickness of the drift field layer in front of the back contact 

for 2 values of mean diffusion length in the drift field layer, 

4 j. 

for an impurity concentration ratio of 10 between p region 
and base region. X m „ => 0.2//m. 
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Just as the light generated current can be expressed as a sum of several' 
current contributions from different regions of the solar cell, so can the 
junction current be expressed as a sum of several current contributions, 
la the most general case, the contributions do not only originate from different 

regions of the diode, but can also be based on various mechanisms. 

Thus , let; 

j d “ j d,D + J d,T + j d,B ; 

where the second subscripts refer to the diffused, transition and base regions, 

respectively. The most basic mechanism leading to dark diode current is the 

diffusion current, described extensively by Shockley P At least in the diffused 

region, it is necessary to take account of the surface recombination velocity. 

2 ) 

Lindmeyer and Wrigley have discussed the case of the narrow base diode with 
an ohmic contact coveilng the entire surface parallel to the pn junction, that is, 
a surface of practically infinite recombination velocity. The solar cell is more 
complicated; the width of the diffused region is, in current cplls, not small com- 
pared to the diffusion length, although in ultimate good cells this should be 
the case. At present, the surface recombination velocity is high everywhere, 
but in ultimate good cells, this should be the case only under the ohmic 
contact strip and the grid lines . 


Thus, a more detailed theory is needed, but based on the same principles as 

used by the above cited references . Of concern is the continuity equation for 

excess minority carriers: $~p - p with p assumed to be constant: 

o o 


A. A 

££= _ £ + D = 0 . 

dt nr p Q ^ ' (US" 

p is the total, instantaneous minority carrier concentration, p^ that which is 
established in thermal equilibrium conditions . p may be larger or smaller than 
P q , with excess carriers recombining or being generated, respectively, with 
time constant n? • The continuity equation ts equal to zero for the steady state 
case, to be considered here. The general solution for the differential equation 


A 

P = Pj e 


+ P 2 e 


with L = L/D 
P P 
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(III-2) 






Using the boundary conditions lor the diffused, layer: 

\A1 

A 


■ d p a 
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= -s p (0) ; 
x =o 


-P (xj = p : 


D 

2_ 
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X' 
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Po — 


B -SPi 


X, 


_ L 

Pj_e P + P 2 e 


p = p ; 


“■ V 


This leads to the matrix: 
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P 

x , 


+s 


L n 
e P 



and the solution of the main determinant: 

x. 
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(III— 3a) 
(III- 3b) 
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(ni-s) 
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--LcoshJL- + s sinh-£— 

X p L P L P . 

D x\ " ^7 

2_ cosh__i + s sinh 

t T, !*„ 


The current into the transition region is thus: 


d,D = -q D 


x=x. 

) 


D 

-q — P- P 


D n x. x. 

^sinh^J + s cosh —L 
L P L P L p 

D x. x. 

— P cosh — - + s sinh - — 
L p L P L P 


° r: ^d, D ~ r 




' <2Y_ ) 

ekT -1 h' 


where; F= 


, , X- sLt, , x i 

sinh , L + - P cosh . ... 

L P L p 

X; sL n Xj 

cosiuJ_ + __JL sinh — — 

^P D p L p 


The factor F includes the dependence on layer thickness^ while the remainder 
of the equation is the familiar relationship for the dark diode current from in- 
finitely thick regions. In fact, the factor F becomes 1 for " 

For SL P^>1, the F becomes cosh-?- , yielding the expression of Lindmeyer 

^ * P 

and Wrigley for the narrow base diode. Fo r P <^1, on tbe °'d 1 ' er band, the 

x,. D p xj 

factor F becomes tanh-£^ , which is less than 1 for small values of*^ . 

Thus with thin regions and low surface recombination velocity , the diode 
(or junction) current can actually be reduced below the normally used values. 
Since reduction of surface recombination velocity has already been found to be 
a goal from the collection efficiency viewpoint, as well as making , it is 

satisfying to see that the same requirements are also beneficial for improvement 
of the IV characteristic. 


- 41 - 



Detail data, as they will be needed for design purposes, of the dependance 
of this factor F on layer thickness and the ratio of surface recombination vel- 
ocity to ^ are given in Fig .HI ~-'l .and UI - 2 . 


It is clear that the ohmic contact areas , including the grid lines , have 
high effective surface recombination . In present cells which have high sur- 
face recombination everywhere, this is not of much significance. In cells 
with reduced surface recombination velocity and low values of ^ > as 

needed for high collection effeciency in the short wavelength region, the 
contact and grid line areas could contribute greatly to the junction current. 

To avoid. this it will become necessary to diffuse deeply below the contact 
and grid line areas, as illustrated in Fig. HI -Since it will (generally not be 
desirable to see a contribution to the junction Current from the contact and 
grid line areas greater than the contribution from the remaining cell surface 
this will provide a minimum val^ for the junction depth x. # c under the 
contact areas. Designers may actually want to use lower contributions 
than this . The total solar ^Tfjuncttori current contribution from the dif- 
fused layer will then be: 



where Ag is the area of the open surface of the diffused region and 
the area covered by contacts and grid lines , including any edge or "wrap 
around" areas , and where F and F "are the F factors applying to A g and A c , 

respectively. p Q r the designer's convenience, a curve of F^, versus _L I s included 

in Fig . Ill - 1* labelled sL p = ^ # 

D 

P 
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Contact 


Diffused Region 




Example: 

x i 

If for collection efficiency reasons, — ^ = 0.1 is chosen, with 

L 

2 P 

s =1000 cm/s. Dp = 1~- and L p = 3.10 4 cm (?= 9-lff 8 s), then F g 

will be 0.38. If is 10% of the total junction area, then Fq should not 

exceed 3.4, requiring x. ^0.3 L p under the contacts. For s = 100 cm s” 1 , 

the corresponding situation leads to^j_ C 

L . 

P P 

It is thus seen that the deep diffusion requirement is not very severe for many 
cases Which prove adequate from the collection efficiency viewpoint. It is 
recognized tnat it will be necessary to extend the deep diffusion sideways by 
an equal distance x^ ^ in order to realize the full effect on the junction current 
contributions . 

This deep diffused region, not covered by contact metal and thus exposed to 
incident light radiation, has a lower collection efficiency than the shallow dif- 
fused regions. Including grid lines, this exposed deep diffused area may comprise 
as much as 5% of the exposed cell area, resulting in a 1% decrease in light generated 
current from the cell if the collection efficiency is reduced by 20% in the deep diffused 
regions. - * 

Since the single layer model for the dark diode current treated so far is not really 
representative of the 2 - layer model, including drift fields-, -used for the collection 
efficiency calculations, a corresponding 2 - layer model lias been developed for the dark 
diode current (saturation current). First it was attempted to reduce the collection ef- 
ficiency equation ( equation ( 33 ) of ref .3) to the zero light intensity case, but this 
was not found to be readily and uniquely possible. In consequence, a new equation 
was derived in complete analogy to the derivation of saturation current in the first part 
of this section and to that of collection efficiency m ref .3 . The new equation was ex- 
pressed in WATFIV language for ease of evaluation by computer. 

Subsequently, a multivariable experiment was run on the computer. Although time 
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did not permit complete evaluation of the data, it appears that, in general, measures 
that improve collection efficiency do also reduce the saturation current, and vice versa. 

A further description of the derivation and an analysis of the results will be in- 
cluded in the next semi annual report. 
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IV 


EVALUATION OF EXPERIMENTAL 


CURRENT-VOLTAGE DATA 


1. Iritroduction 

To carry out optimum-design studies relating to the current- voltage 
characteristic, it is imperative that the physical theory determining this 
characteristic be understood. This understanding has, at least, to be 
sufficient to name all physical and material parameters which do actually 
have influence on the IV characteristic of silicon solar cells, and to 
express the relationship of these parameters to the IV- characteristic. 

This understanding does not yet exist, although there are a number of 
theories on pn-junction behaviour. What is first required is accurate 
measurement of current- voltage characteristics under various conditions 
on numerous solar cells, and comparison of the experimental data with 
the prediction of theory. This section reports on some further work to- 
wards this aim. 

As discussed in section IV of the previous semi-annual report, the 
current-voltage characteristic of silicon solar cells is better described 
by the relationship 

j = j 0 l - 1 1 + J 0 2 [ eXp ^A 2 kT^ “ 1 1 ; (IV_1) 

( 3 is used for current density, all other symbols as used before ) 
than by the more commonly applied one containing only one such expon- 
ential term. Because of its mathematical features, this relationship will 
from here on, be called the "two-exponential IV-characteristic" , in con- 
trast to the common form, which may be referred to as the "single expon- 
ential IV-characteristic". The following subsection discusses a number 
of observations made in the interpretation of experimental data of IV 
characteristics, and provides some useful hints to other workers in this field. 
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2 . On the interpretation of experimental data of XV- characteristlcs^ 

1 

A few interesting observations have been made during the course of 
this author's work on the current- voltage characteristics, which may be 
of considerable benefit to the readers if shared with them. These observa- 
tions are concerned with the use and interpretation of the experimental data 
on the current- voltage characteristics of solar cells , and are only in secon- 
dary manner connected with the computer program for their evaluation. 

Figure IV-1 shows a set of experimental current voltage points (solid 
dots) in the range above 0.25Volt. The classical approach in the interpre- 
tation of such data has been to draw two tangents through the straight line 
portions of the characteristic, and by determining the slope of these tangents, 
compute the value of the factors and A 2 in the exponents of the two-exponen- 
tial IV characteristics. This author has, for some time, found that this approach 
generally does not lead to satisfactory results. The superposition of the 
two exponentials with the thus found A values does not usually represent 
the experimental points , particularly where they fall into the overlap 
region between dominance of one of the two exponentials, or into the 
region where the -I Q2 term has significant influence. To overcome this 
shortcoming, the author has used a manual iterative method which is illus- 
trated in Fig. IV-1. Here, two straight lines are shown as tangents to the 
experimental data, labelled A 2 =8 .46 andA^l.24. The exponential term with 
the latter A value fits to the experimental data in the region above 0.55Volt, 
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with a series resistance value of .04 Ohm. In the overlap region between 
the two exponentials, substracting the values of this exponential term from 
the experimental data points should lead to values represented by the second 
exponential term containing . The resulting points are indicated by the 
crosses in Fig. IV- 1. It is clearly seen that above 0.4Volt, the two exponen- 
tial terms discussed so far do not represent the experimental data, since 
the crosses fall off towards smaller current values at higher voltage, rather 
than continuing on the straight line representing the exponential term with 
A 2 =8.46. A change in the slope of the exponential with A^ can remedy 
this situation as indicated by a second straight line, labelled A »1. 12 in 
the semi-logarithmic plot. Subtracting the values of this exponential term 
from the experimental data points leads to new current values in the region 
between 0.4 and O.SSVolt. As seen, these values fall nicely onto the 
straight line of the exponential term with A 2 =8.46. Thus, the experimental 
values can.be represented by the superposition of two exponential terms 
with A =1.12 and A =8. 4 6, but not by use of an exponential term with 
A^=l . 24 . The latter, however, would have formed a tangent to the upper 
part of the experimental curve , while the exponential term with the smaller 
A^value intersects the straight line portion of the experimental curve rather 
than forming a tangent. This observation has been made repeatedly, and 
will be substantiated with further detail in the following paragraphs. 

Figure IV-2 shows the same set of experimental data, this time however, 
extended to voltages below O.lVolt. It is seen that the lowest two experi- 
mental points fall onto the straight line with slope A 2 =8.45. However, 
an effect of the -I term would have been expected to be noticeable on 

\jc* 

these lower points. The extrapolation of the exponential term to its inter- 

-5 

sect with the ordinate yields the value of I^ 2 , Using the value of 7*10 Acm 
obtained in Fig. IV-2, and subtracting this value from the experimental data 
yields the curve represented by crosses. It is seen that this curve differs 






considerably from the actual experimental data. It is also seen that any 
attempt to use an exponential term with a smaller slope to explain the 
experimental data in the region above 0.2Volt, will lead to a larger value 
of Iq 2 and, consequently, to a larger curved part of the resulting charac- 
teristic at the lower voltage values. As a compromise, an exponential 
close to that labelled A 2 =12.47 has been chosen, and the resulting satur- 
ation current I has been subtracted from the points on the straight line 
representing this exponential term in the semilogarithmic plot. The result 
was a significant deviation from the experimental data only at the two 

lowest points, as indicated by the hollow circles. 

It is interesting to note that the computer program VICHAR was not able 

to converge with the original experimental data points. It was able, however, 
to find convergence for the curve composed of the experimental data points 
(solid dots) above 0.46Volt and of the corrected points marked by crosses, 
at lower voltages, as well as for the curve given by the experimental data 

points above 0.2Volt and the hollow circles below. 

It was later determined that the experimental data had been faulty. 

It is very easy to obtain erroneous data points from the instrumentation 
at such low current values due to inappropriate internal resistance of the 
measuring equipment, due to noise pick-up, or ground loops. Thus, inability 
to fit all data points to a theoretical relationship, or large values of standard 
deviation obtained by use of the computer program VICHAR , or inability 
to converge may well be caused by faulty experimental data , and should 
lead to a re-examination of the instrumentation if such occurrances are noted. 

In order to further analyze the relationship between the values of the 
individual exponential terms and the points resulting from their superposition 
with the inclusion of series resistance, three models of solar cell character- 
istics were prepared. A second purpose of this undertaking was to check 
the ability of the computer program to find the actual values which were 
used to generate such a curve. For Model 1, the values of A =1.0, A 9 =2.0, 

0 o 0 2 x 

j0-|=lO A cm , and 3 o2 == '*‘ ^ A cm were chosen, together, with R=0.5 Ohm, 
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MODEL 1 

Two- Exponential Current-Voltage Characteristic 
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^1= 1.0; 

I Q2 = 1.10 6 Acm 

■ 2 ; A 2 = 2.0; 

R = 0.5 

yi 


V D 

I Q1 e xp ( 

qV D ) 

I n o ex P ( 

qV D ) 

— ^fi? 

I 

\ 

V 

A i kT 


A 2 kT 





0.01792 



1.4142 

•io -6 

- 1.0 *io“ 6 

4.142 

• io -7 


0.03584 



2.0 


-1.0 

1.0 

•io- 6 


0.08321 



5.0 


-1.0 

4.0 



0.11904 



1.0 

•10" 5 

-1.0 

9.0 



0.15488 

1.0 

• io“ 8 

2.0 


-1.0 

1.901 

• io- 5 


0.2023 

2.5 

• 10" 7 

5.0 


-1.0 

4.925 



0.2381 

1.0 

• 10- 6 

1.0 

• io' 4 

-1.0 

1.000 

♦ io -4 


0.2739 

4.0 


2.0 


-1.0 

2.03 


0.2740 

0.3213 

2.5 

• 10 -5 

5.0 


-1.0 

5.24 


0.3216 

0.3571 

1.0 

• io" 4 

1.0 

• io -3 

-1.0 

1.099 

* io" 3 

0.3576 

0.3930 

4.0 


2.0 


-1.0 

2.399 


0.3942 

0.4167 

1.0 

• 10 -3 

3.162 


-1.0 

4.161 


0.4188 

0.4346 

2.0 


4.472 


-1.0 

6.471 


0.4378 

0.4583 

5.0 


7.071 

• 

-1.0 

1.207 


0.4643 

0.4762 

1.0 

• io -2 

1.0 

* io -2 


2.0 


0.4862 

0.4941 

2.0 


1.414 



3.414 


0.5112 

0.5178 

5.0 


2.236 



7.236 


0.5540 

0.5357 

1.0 

• io -1 

3.162 


• 

1.316 

* io -1 

0.6015 

0.5536 

2.0 


4.472 



2.447 


0.6760 

0.5773 

5.0 


7.071 



5.707 


0.8627 


TABLE IV-1 
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T=300°K. The data for current density were computed to four digit accurracy, 
and the voltage values were corrected, where appropriate, to reflect the 
voltage drop across the series resistance, to the same accurracy. The values 
of the exponential terms themselves, the j Q2 value, their superposition, 
and the corrected voltage values are shown in Table IV- 1 . The resulting 
current- voltage points are plotted as solid dots in Fig. IV- 3. The result- 
ing IV characteristic seems to have straight line portions over more than 
two orders of magnitude in the lower voltage region and of approximately one 
order of magnitude in the higher voltage region. The tangents to the straight 
line portions of the resulting IV characteristics are shown as dashed lines. 
However, the A values for these straight line portions are A 2 =1.87 and 
A^=l . 67 , while the values generating the data were A^l . 0 and A 2 =2 . 0 . 

These original exponential terms are represented by solid lines in Fig. IV-3. 

It is seen, that the exponential term containing A^ is far from being tangent 
to the straight line portion in the upper part of the characteristic, but, 
rather intersects it at an angle of approximately 30°. It is also clear that 
quite a bit of manual iteration would be required to find the appropriate expon- 
ential terms to explain the experimental data. This author has frequently 

spent two or three hours for the manual iteration process on one set of 
' 1 
experimental data. 

Table IV-2 shows another set of data, called Model 2. The only differ- 

-7 -2 

ence in comparison to Model 1 is a change to j q 2 = 1° A. cm . Figure IV- 4 

presents the corresponding data points, the tangents to the straight line 

portions of the characteristics, shown by dashed lines having A^=1.19-, and 

A =1.84, as well as the generating curves with A =1.0 and A =2.0. While 
Z 1 z 

the table shows the large overlap regions between the various terms in the 
equation, the figure indicates straight line portions of the characteristic of 
approximately two orders of magnitude for both the A^ and the A 2 regions . 
Again, the generating exponential terms are represented by intersecting lines 
rather than tangents. 
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Current Density 


Solar Cell Model 1 

j ol = l.lO" ,0 A cm" 2 ; A, = 1.0; 
j 02 =I.IO' 6 A cm" 2 ; A 2 = 2.0; 
R = 0.5X1 ; T = 300°K 
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Voltage V 


The current-voltage characteristic of Model-1 , together with 
the two straight lines representing its generating terms (A =1.0 
and * 2 = 2 . 0 ), as well as with tangents to the straight line portions 
of the characteristic. 





MODEL 2 


Two-Exponential Current-Voltage Characteristic 
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Voltage V 


Fig. IV- 4 The current-voltage characteristic of Model-2, together with 

the two straight lines representing its generating terms 
(A = 1.0 and A 2 = 2.0), as well as with tangents to the straight 
lino portions of the characteristic. 





Subsection 3 then reports further work on a computer program (called 
VI CHAR for voltage- current characteristic) to ease the interpretation 
of such experimental data. The mathematical background for this 
program was provided in section IV of the previous semi-annual report. 
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It is thus evident that considerable care needs to be taken, both in 
the acquisition of the experimental data and in their interpretation. Manual 
interpretation of the two-exponential characteristic leads to a time con- 
suming iteration process. Here, the computer program VI CHAR has proven 
to be an extremely valuable tool , providing not only considerable time 
savings, but also an accuracy in data evaluation not matchable in the 
manual process. More on this subject is contained in the next subsection. 
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3. The Computer Program VIGHAR 

In section IV of the first semi-annual report, work-in-progress on a 
computer program for the evaluation of the five constants , (including 
series resistance) in the "two-exponential" current voltage characteristics 
of solar cells has been discussed. Further work on this program (called 
VI CHAR) has been performed in the present reporting period, with the 
result that considerable confidence in the accuracy and usefulness of 
this approach have been gained. As a consequence / the- work ©n the 
generation of this program can be considered completed, with its extensive 
application as a tool in the future research on current- voltage character- 
istics indicated. 

A program of the mentioned purpose is really useful only, if it will 
arrive at correct answers for every realistic case of experimental data, 
without any actions by the computer operator beyond initiating the input. 
This condition was not fulfilled for a long time due to the fact that the 
approach is rather sensitive to deviations from their true values of the 
constants in the exponents. This is particularly the case for R, the 
value of the series resistance. As long as the estimated (or starting) 
values for the five constants , or their corrected values in the early iteration 
steps, are not reasonably close to the ultimate values represented by the 
experimental data, there is a tendency in the program to over-correct, 
resulting in possible divergence rather than convergence, to arrive at 
exponents greater than provided for in the computer , (particularly at low 
temperatures ), or at negative values of R or the j q* s . A substantial number 
of safeguards have therefore now been built into the program, which has 
resulted in considerably greater complexity than contained in the original 
version. 
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As one of these safeguards , the program hunts for the most rapid 

convergence by cutting all the correction values resulting from solution 

of the matrix sequentially in half, down to 2 Similarly, it avoids 

using negative values of resistance or saturation current. Further, the 

program limits the exponents in magnitude to that permitted in the IBM 
+7 fi 

360 system (10 ), and if it cannot find convergence with an estimated 

starting value for series resistance, it will try with up to 5 other values. 

As a result, the program has lately run to completion on all sets of input 
data, which reasonably closely represent a two- exponential diode charac- 
teristic as expressed by equations 2,5, and 6 in section IV of the first 
semi-annual report. The program has proven rather sensitive to this 
condition, rejecting data containing erroneous points, but has ’found 
convergence on input data containing only one exponential or R=0. 

Figure IV-5 presents an abridged flow chart of the present program, 
showing its essential features. VTCHAR can be used to determine either 
all 5 constants, 4 constants, by entering predetermined values for either 
series resistance R, as may be obtained from separate measurement, or 
for A^ ,( called Al in program language (WATFIV)),or 3 constants, by enter- 
ing values for both R and A^. Accordingly, the program has to solve systems 
of 5, 4, or 3 simultaneous equations, respectively. For this purpose, 
either 30 constants ( 5 by 6 matrix ), 20 constants ( 4 by 5 matrix ) or 
12 constants will be computed. The 3 to 5 correction values ORANS (y) are 
found by solution of the simultaneous equations through diagonalization 
of the matrix in subroutine SIMUL . The quality of fit between the experi- 
mental data and the theoretical two-exponential characteristic at any stage 
in the iteration process is determined by the second moment of deviation 
(sum of the squares), evaluated in function ERR, Application of the correction 

values to either the estimated (or starting) values (ending with - EST) for 
the 3 to 5 constants to be determined, or to the values of the constants 
resulting from prior corrections, is made in subroutine USANS. 
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Most of the remainder of the program is devoted to assuring convergence, 
recycling, etc. , as discussed previously. Finally, the program in its 
present form , after having been run to determine a value for , will re- 
cycle with the fixed value A^= 1. 0. 

The convergence criteria for ending the iterations are: no relative 
change in any of the constants by more than 0.0005 in the last iteration 
step, or, after 30 iteration steps, the standard deviation between all 
input data and the two- exponential characteristic with the computed 
constants is less than 10% and has not changed more than 0.0001 in 
the last iteration step. 

It may be noted that all effort in developing the program has been 
directed towards its performing the required function, without any 
attention paid to an elegant approach or lowest cost execution. 


As a significant test of the program VI CHAR , the IV data of Model 1 
and Model 2, discussed in subsection |V- 2 above, have been subjected 


A 2 was 1.9995, with j Q2 = 


to its evaluation. On Model 1, the program took 22 iteration steps to 
arrive at A =0.9993 and j =0.985*10 

X V X 

0.9989 • 10 , and R=0.5001. The standard deviation between the entered data 

-4 

points and the theoretical curve having the resulting parameters , was 7*10 
With the exception of j , all values are accurate to the fourth significant digit. 
Similarly, on Model 2, the data were reproduced to the fourth significant 

digit. Again, j had a slightly larger error, and the standard deviation was 

-4 U1 

8*10 in this case. 


Two other tests were run on the models, one of them using Model 1 
data without series resistance. This is equivalent to using junction data 
instead of terminal data. The program came to convergence within 10 
iteration steps , again finding the original parameters to the fourth signifi- 
cant figure except for j , which had a 2% error ,and with a standard de- 
-4 U1 

viation of 5*10 . In the other test, only one exponential term, the one 
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with A =1.0, was used, but including series resistance R=0.5 Ohm, (called 

Model 3). Here, the program surpressed j Q2 to a sufficiently small value, 

to eliminate the influence of the second exponential term. Again, the 

values of , A, , and R were reproduced to better than four significant 
01 1 

digits , with a standard deviation of 9 * 1 0 

These tests withdhe model data have shown clearly that the program 
VI CHAR is capable of reproducing the five constants which were used 
to generate the input data, with essentially the accuracy of the original 
data. Thus, the program does not only provide a speedy tool for the evalu- 
ation of the five constants, but, moreover, provides an accuracy which 
cannot be matched by the manual iteration process. 

An additional interesting test was accidentally performed, which 

has shed more light on the question of interpretation of erroneous data. 

In preparing the card deck for Model 1 without series resistance, one 

card from 'the deck with series resistance was left in the new deck. It 

-5 

had a relative error in voltage of approximately 4.10 compared to the 

newcard, and was left out of sequence in the deck. Both of these facts 

should not influence the results < significantly r- Further, a keypunch error 

had occurred on the 34.14 m A current value, resulting in 31.0m A. 

The error was corrected, but the faulty card left in the deck. Since the 

deck had 2 cards more than expected, the computer did not use -the last 

2 cards, for current values of 244.7 and 510.7 m A. As a result, A^ 

came out 5.7% too high, and j was high by a factor of 2.53. A 2 and j Q2 

differed only in the 4th significant digit, while R was correct. The standard 

-4 

deviation was 4.4%, instead of 4 . 7* 10 obtained with the correct deck. 

Instructions for use of the program VT CHAR are included as Appendix I 
to this report. Figures IV- 6 and IV-7 show two typical pages of computer 
printout, the first listing the instructions given and the input data, as well 
as their correction to dark diode junction data , using either an entered value 


- 66 - 



CEVICE NUMBER =' MCDEL 1 
AREA NET GIVEN 

PRINTOUTS IN A FOR I, OHM FOR R 

DARK DIODE, TERMINAL DATA 
RESISTANCE UNKNOWN 

TEMPERATURE = 300. CO 

NUMBER OF DATA PAIRS = 20 


VK IK 


C.17920D- 

•Cl 

0.41420C-C6 

C.35840D- 

•01 

0 . 1CG0GD-C5 

C.83210D- 

■Cl 

0.4CCC0D-C5 

G. 1 19040 

CO 

0.9CCC0D-C5 

C.15488D 

CO 

0. 19010D-C4 

0.2C230D 

00 

0.492500-C4 

0 .238100 

cc 

0. 10C00D-C3 

0.274000 

cc 

0.2C3CC0-C3 

0.32160D 

cc 

0.524G0D-C3 

0.357600 

CO 

0 • 1C990D-C2 

0.3 94200 

cc 

0 . 239900- C2 

0.41880D 

CO 

0 .41610D-C2 

0 . A 37800 

CO 

0.6472GD-C2 

C.46430D 

cc 

0. 12070D-C1 

0.466200 

CO 

0.20CGCD-C1 

C.51120D 

cc 

0 .34140D-C1 

0.55400D 

cc 

0 . 723 600- C 1 

C.6C150D 

cc 

0.13160D CC 

0.676000 

cc 

0 .244700 CC 

C.86270D 

cc 

G . 570700 CC 


VDK 


I OK 

0.179200- 

01 

0.414200-06 

0.35839D- 

01 

0. 100000-05 

0. 832C6D- 

01 

C . 4CGOCD-05 

0.119030 

CC 

0 .900000-05 

0.154860 

00 

0. 190100-04 

0.202250 

CC 

C. 492500-04 

0.238000 

OG 

C. 1000CC-C3 

0.273800 

CO 

C. 203000-03 

0.321C8D 

00 

0.524000-03 

0.356500 

CC 

0. 1099CD-02 

0.391800 

00 

0.239900-02 

0.414640 

00 

0. 416100-02 

0.431330 

CC 

0. 6472CC-02 

0. 452230 

OG 

C • 12C7CD-01 

0.466200 

CO 

0.200000-01 

0.477060 

CC 

C. 34 14CD-G1 

0.48 1640 

CO 

G . 7236CD-01 

0.469900 

00 

0 . 13 1600 CO 

0.43130D 

OG 

0.244700 00 

0.292C0D 

CO 

G.57070D CC 


IV- 6 Printout page resulting from program VICHAR, showing the input 
data, in this case for Model -1. 
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ITER »N 

IC1 

A 1 

1D2 

A 2 

R 

SIGNA 


IMT'L 

0.5CCCCD-1C 

1.CCCC0 

C • 2CG00D-05 

2.5CC00 

l.COOOO 

0.6C946D 

OC 

1 

0.645 9 5D- 1 C 

1.01722 

C. 199130-05 

2.49654 

0.96415 

0.60768D 

OC 

2 

0. 8279 4 D- 1C 

1.03431 

C . 198280-05 

2.49317 

0.92962 

0.60573D 

OC 

3 

0. 1C5220-C9 

1.05115 

C . 19745D-05 

2.48988 

0.89655 

0.6C3600 

OC 

4 

0. 13248D-C9 

1.06766 

C . 196620-05 

2.48668 

0.86505 

0.60130D 

OC 

5 

0. 19787D-C9 

1. 1C029 

C. 19500D-05 

2.48044 

0.80538 

0.59859D 

00 

6 

0.92815D-C9 

1.41787 

C . 182380-05 

2.43355 

0.36915 

0.59732D 

OC 

7 

0.92815D-C9 

1.19540 

C. 180740-05 

2.42230 

0.41979 

0.420800 

00 

8 

0.128150-C8 

1.15656 

-C. 13169D-05 

2.21156 

0.47765 

0.16710D 

00 

9 

0.97843D-C9 

1.13928 

C . 12693D-05 

2.17999 

0.48055 

0.156560 

OC 

1C 

0.7549 1D-C9 

1.12302 

C. 122950-05 

2.15356 

0.48314 

0 . 1467 ID 

00 

11 

0.59213D-C9 

1.10801 

C. 119620-05 

2.13138 

0.48544 

0.13702D 

OC 

12 

0.47-401D-C9 

1.09439 

C . 116820-05 

2.11268 

0.48745 

0.127270 

OC 

13 

0 . 388CCD-C9 

1.08220 

C .114450-05 

2.09687 

0.48920 

0. 11742D 

OC 

14 

C.3248CD-C9 

1.07141 

C. 1 1 244D-C5 

2.08344 

0.49071 

0.107600 

00 

15 

0.27 78CO-C9 

1.06194 

C .110730-05 

2.07201 

0.49202 

0.979660- 

-01 

16 

0.20692U-C9_ 

1.04556 

C. 107810-05 

2.05255 

0.49425 

0.867260- 

-01 

17 

0.1 6748D-C9 

1.03314 

C. 105710-05 

2.03856 

0.4.9588 

C. 721250- 

-01 

18 

0. 12125D-C9 

1.01514 

C. 10267D-05 

2.01833 

0.49821 

0 . 52619D- 

-01 

19 

0.94226D-1C 

0.99825 

C .99829D-06 

1.99924 

0.50027 

0.988310- 

-02 

20 

O'. 984870- 1C 

0.99934 

C .898880-06 

1.99946 

0.5C009 

0.762240- 

-03 

21 

0 . 98 507D- 1 C 

0.99931 

C .998880-06 

1.99946 

0.50010 

0. ,704710- 

-03 

22 

0 .985070- 1C 

0.99931 

C. 898880-06 

1.99946 

0. ; 50010 

0.7C471D- 

-03 


IV- 7 Printout page resulting from program VICHAR, showing the results, 
in this case for Model-1. 
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for series resistance R or the estimated value REST, and an entered value 
for the light generated current IL, if applicable. Fig. IV- 7 contains the 
results of each iteration step, starting with the estimated values in the 
first line, and ending with the final values for the 5 constants, after the 
final convergence criteria have been satisfied. The first column contains 
the number of the iteration steps before printout of that line, and the 
last column shows the standard deviation SIGMA achieved with the con- 
stants shown on the same line in the remaining columns. 

The discussions relative to Fig. IV-2 have shown that 2 data points 
off the characteristic can cause no-convergence. The accidental test 
with the faulty deck, discussed above, showed that a 10% error in one 
data point can cause an increase of about 4% in the standard 
deviation. These are indications of the great sensitivity of the method 
to errors in the input. From other tests, the conclusion has been reached 
that results with standard deviations above 10% should definitely be con- 
sidered as insufficient convergence, requiring examination of the input 

data. Also, the input data shall span from at least 50mV to above 
-2 

lOOmAcm , with at least 2, but preferably '3 data. points per decade. 
Maximum errors in input data of 0.1% for voltage and 1% for current should 
be strived for. It may be noted that a 0.1% error in voltage at 0.6V 
corresponds to a 2.4% current error. Combined with a 1% current error 
(independent measurement), this means a maximum deviation from the 
true value of 3 . 4% , or a probable deviation of 1.7%. Thus , a standard 
deviation for the entire data group around 1.5% might be expected. 
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4. Conclusion: 


The program VICHAR has proven to be an excellent tool for the 
evaluation of experimental data of the current- voltage characteristic ' 
of silicon solar cells, and possibly other pn-junction devices. It 
does not only provide a tremendous labor saving compared to the 
earlier manual method, but, also yields much- greater .accuracy-.} 
provided the input data are of sufficient .quality. As the result of the 
availability of this new tool for the evaluation of the IV characteristics, 
it has become clear that further studies will require more accurate experi- 
mental data, than have generally been used up to now . However, one 
conclusion can already be drawn from the limited amount of experimental 
data evaluated with the computer program: the small standard deviations 
(2 to 4%), with which the two-exponential-* current- voltage relationship 
has been fitted to the experimental curves, strongly suggests that the 
two-exponential characteristic actually is the analytical relationship 
which describes the experimental data. Up to this point, the doubt 
existed that the two-exponential characteristic might be nothing but 
a reasonably good approximation to the experimental data, which might 
otherwise be better described by a different mathematical relationship. 
However, with this close reproduction of the experimental data, and 
in view of the great sensitivity to data deviations discussed earlier# 
considerable confidence exists -now that the -two-exponential relation- 
ship is actually the proper analytical expression for the IV character- 
istic of silicon solar cells and perhaps for many other pn-junction devices. 

In consequence, the next step will be to re-examine the existing 
theories on pn-junction current flow for their capability to provide the 
required two-exponential IV characteristic. If they should not be capable 
of describing this characteristic, the search would have to go much 
deeper, into basic phenomena and their relationship to the IV- characteristic. 


- 70 - 



V IMPROVEMENT OF THE DIFFUSED REGION 

COLLECTION EFFICIENCY 

* p* 

It had been predicted (ref.l), that a 17% improvement in collection 
efficiency should be obtainable by reduction of the surface recombination 
velocity of the diffused region and by reduction of the diffused region 
thickness. Later design calculations, included in the first semi-annual 
report, indicated that the diffused region thickness reduction might not 
be necessary, if the quality of the material parameters in the diffused 
region should be better than originally thought, or if it could be improved 
by a relatively small margin. It was expected, that some preliminary 
experimental work towards the verification of this prediction would be 
carried out. 

As the planning for this experimental program got under way at the 
beginning of this year, it was learned accidentally, that. a program with 
the same objective, spawned by this author's prediction, was well 
under way at Comsat Laboratories uner the direction of Dr. J. Lindmayer. 

It was also learned that preliminary results appeared hopeful. Since 
more extensive and conclusive results were expected within a few months, 
the preparations for our own experimentation were slowed down, in order 
to avoid unnecessary duplication. 

It has now been learned from Dr. Lindmayer, that the experiments 
have been very successful. The collection efficiency of the silicon cells 
has been improved by about 20%, based, however, on a cell of thinner 
base region thickness than used in the calculations. The lower base 
region thickness implies a lower long wavelength collection efficiency 
and thus, a lower total collection efficiency base for the improvement. 
Further, the anti-reflection coating was heavily blue shifted, with a 
resulting higher reflection in the long wavelength region, resulting in 
an additional effect of the same consequences as the reduced base re- 
gion thickness. Thus, the achieved improvement is just of the right 
magnitude . 
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Details of the method by which the improvement was achieved are 
still proprietory information of Comsat Laboratories, but are in prepara- 
tion for publication. It has also been found there, that the conventional 
process providing the anti-reflection coating does not permit full realiza- 
tion of the improved short wavelength response achieved by the diffused 
region parameter improvement. Further development of suitable anti- 
reflection coatings is still under way at Comsat Laboratories , which has 
a publication in preparation on the subject of the short wavelength response 
improvement. 

It may also be noted with interest, that the collection efficiency 
improvement in the short wavelength region has been found to be ex- 
tremely radiation resistant in preliminary tests at Comsat Laboratories. 
Thus, the improvement in overall collection efficiency of the silicon solar 
cell has been accompanied by a new level of radiation hardness. This also 
is in line with the expectations. 

In view of this work performed at Comsat Laboratories, it would not 
appear sensible to proceed with the program as originally planned. Rather, 
further disclosure of data from Comsat Laboratories should be awaited. 
Following evaluation of these data , a modified program of experimental 
investigation may then be undertaken. 
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APPENDIX 1 


Instructions for Use of Program VI CHAR 


1 , Cards of Data Package: 

Start card: /GO in columns 1 to 3 
Card 1: Estimated or starting values for 

I 01 (IjariEST ), A X (A1EST), I 02 ( IjefEST ), A 2 (A2EST), 

_2 

and R ( REST ), in this sequence. I in Acm or A, R in Ohm. 
FORMAT : 

Iq ^ = double precision floating decimal ( D20.14 ), 
start column 1; 

A^ = fixed decimal ( F10.6 ), start column 21; 

I = double precision floating decimal ( D20.14 ), 

Ua 

start column 31; 

.A^ = fixed decimal ( FI 0.6 ), start column 51; 

R = fixed decimal, (FI 0.6 ), start column 61. 

Card 2: Number (NSYS) of VI - characteristics to be processed. 

FORMAT: Integer, columns 1-5 (units in column 5, 
tens in column 4, etc.) 

The following cards have to be provided for 
for each VI - characteristic to be processed. 

Card 3: Device temperature T, in degrees K; 

FORMAT: fixed decimal, columns 1-20 

Card 4: Number of fixed resistance values to be used (NRES), fixed 

resistance values (RES) in Ohm; (up to 15 values). 

FORMAT : NRES: integer columns 4-5 (units in column 5, 
tens in 4, etc.) 

RES: fixed decimal (F10.6), 7 times, starting column 11, 

additional 8 times on continuation card, starting 



Card 5: Light generated current I (IL) , in Acm or A, when 

L 

illuminated VI- character! Stic is to be evaluated. Otherwise 
no entry on card. 

FORMAT: fixed decimal , columns 1-20. 

Card 6: Device label (identification) (DEV), number of data pairs 

2 

in Vi-characteristic (ND), junction area in cm (AREA), 
in this sequence. 

FORMAT: alphameric columns 1 to 16 for DEV; integer 

columns 21 and 22 for ND (units in column 22); 
double precision floating decimal for AREA, columns 
31 to 50. 

-2 

ND cards: voltage (V) in Volt and current (I) in Acm or A, in this sequence, 
FORMAT: V: fixed decimal columns 1 to 20, 

I: fixed decimal columns 21-40. 

Each data pair one card. 

End Card: / END columns 1 to 4. 

_2 

2. All current values to be entered either in Acm or in A. 

3 . A value for AREA should be entered , when current data ( in A) are entered 

-2 

and evaluation in current density (in Acm ) is desired (normalization 
to unit area) . 

4. If terminal voltage data are entered, at least one value for either 
REST or RES has to be entered. Otherwise, the input data will be 
treated as junction voltage data. 

5., If evaluation with a fixed resistance value is to be carried out, no 
value (or 0.0) is to be entered for REST. 

6. If evaluation with a fixed value of 1.0 for A is to executed, no value 
(or 0.0) for AIEST is to be entered. 

7. Any combination of the above is permitted. 
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VI PROTECT STATUS 

1. Point 1 of the program is completed. 

2. ' Point 2 has been delayed as explained in section V of this report. 

3. Point 3 is in progress. Contacts to material suppliers have been 
made and answers are starting to be received. 

4 . As additional work , the computer program VICHAR for evaluation of 
the 5 constants is the 2-exponential current-voltage characteristic 
is completed. 

5 . As additional work , a study of the saturation current contribution 
from the base region, using a 2-layer model with drift fields, similar 

• to the erne evaluated for collection efficiency ,1s well under way. 

6. As additional work, the data resulting from point 1 of the milestone 
chart have been re-cast to form the beginning of a "Silicon Solar Cell 
Design Handbook" . 


- 75 - 



PLANS FOR NEXT REPORTING PERIOD 


Whenever appropriate, redesign experimentation according to 
section point 2 of the milestone chart. 

Proceed with point 3a, of the milestone chart and start point 3b. 

Complete the study of saturation current contribution from the 
base region. 



vm aU.MMA.RYs CONCLUSIONS, RECOMMENDATIONS 


1 . The collection efficiency part of a "Silicon Solar Cell Design Handbook" has 
been completed. For this purpose, the total solar flux is divided into the 
parts absorbed m the various region of the solar cell (diffused region, base 
region, etc .) , (function of region thickness), and the collection efficiency 
from these regions determined as function of the relevant physical and material 
parameters . 

2. Additional data sets for base region collection efficiency were obtained for 
various values of distance from the light exposed surface to the interface plane 
between depletion region and base region. 

3. A relationship was derived to evaluate the saturation current contribution from 
the diffused region for arbitrary surface recombination velocity and arbitrary 
thickness of the diffused region. It was found that a "form-factor" applied to 
the standard saturation current formula best describes the influence of these 
parameters . 

4. Using the relationship mentioned under 3 , it was determined that: 

a. The proximity of the ohmic contact as such (excluding process caused 
effects, such as diffusion of metals) to the depletion region (pn junction) 
has no significant effect in present silicon solar cells due to the high sur 
face recombination velocity at the exposed surface . 

b. In future cells of low surface recombination velocity, saturation current 
below the values predicted from the diffusion theory for thick layers (the 
commonly used relationship) can be obtained. This can result, e.g. , in 
an easing of the requirement for good minority carrier lifetime in the dif- 
fused region. 
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c. If case 4b. , is fulfilled, it may be necessary to provide a deeper diffused 
guard layer under the ohmic contacts . 

5. To investigate the influence on saturation current (and thus on the IV-character- 
istic) of measures found to improve collection efficiency m thin' cells (such as a 
drift field in front of the back contact ), a base layer model for saturation current 
computation was derived in analogy to the 2 - layer model including field terms 
for the collection efficiency. 

6. Using the model of 5 . , above, a multivariable "experiment" was completed, but 
not yet fully evaluated. Preliminary inspection seems to indicate however, 
that whatever improves collection efficiency, also lowers the saturation Current. 

7 . The computer program VICHAR for the evaluation of the 5 constants in the two- 
exponential current -voltage characteristic from experimental data has been com- 
pleted. It not only provides labor saving, but also an accuracy not obtainable by 
other methods . 

8. Experimentation carried out at Comsat Laboratories has achieved the improvement 
of short wavelength collection efficiency predicted from this analytical work. It 
has also been found there that the present anti -ref lection coatings do not permit 
full utilization of this improvement, resulting in efforts there towards development 
of a better AR coating . 
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